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ABSTRACT 



system. The density of states v is expressed formally as a functional integral over 
Lorentzian metrics and is a functional of the geometrical boundary data that are fixed 
in the corresponding action. These boundary data are the thermodynamical extensive 
variables, including the energy and angular momentum of the system. When the 
boundary data are chosen such that the system is described semiclassically by any real 
stationary axisymmetric black hole, then in this same approximation In z/ is shown to 
equal 1/4 the area of the black hole event horizon. The canonical and grand canonical 
partition functions are obtained by integral transforms of v that lead to "imaginary 
time" functional integrals. A general form of the first law of thermodynamics for 
stationary black holes is derived. For the simpler case of nonrelativistic mechanics, 
the density of states is expressed as a real-time functional integral and then used 
to deduce Feynman's imaginary-time functional integral for the canonical partition 
function. 



* Present address: Departments of Physics and Mathematics, North Carolina 
State University, Raleigh, NC 27695-8202 



I. INTRODUCTION 



The energy of a physical system is reflected in the gravitational field it pro- 
duces, so the gravitational field at a (spatially finite or infinite) closed surface that 
bounds the system encodes information about the energy content. By fixing appro- 
priate components of the gravitational field, the energy of a self-gravitating system 
can be specified as boundary data. In statistical mechanics and thermodynamics, 
where the concept of energy plays a central role, this circumstance allows for the 
direct specification of microcanonical boundary conditions in which the thermody- 
namical extensive variables (including energy) are held fixed. Here, we exploit this 
special property of the gravitational field in a direct construction of a "microcanon- 
ical functional integral", a formal functional integral expression for the density of 
states, which characterizes a system with microcanonical boundary conditions. 

The canonical partition function for nonrelativistic mechanics was first ex- 
pressed as an imaginary time functional integral by Feynman.[l] This prescription 
was later generalized to flat space fleld theory, [2] then to self-gravitating systems by 
Gibbons and Hawking. [3] As an alternative to this line of development, we present 
a direct expression of the density of states for nonrelativistic mechanics as a real- 
time functional integral. The generalization of this result to flat space fleld theory 
is not immediate, because in that case flxing the energy involves a restriction on the 
integral of the Hamiltonian over the entire spatial extent of the system. However, 
the generalization to gravitating systems is quite natural, because the presence of 
gravity allows the energy to be flxed directly by the boundary data. In this paper, 
we consider the functional integral expression for the density of states for systems 
consisting only of the gravitational fleld. The inclusion of various matter flelds will be 
given elsewhere. [4] Inasmuch as all systems are self-gravitating, even if only weakly, 
the formalism developed here is in principle completely general. 

One of the key features of the present analysis is the use of flnite boundaries 
in space. There are a number of advantages to be gained by imposing boundary 
conditions at a spatially flnite location, as opposed to spatial inflnity. For example, 
with flnite spatial boundaries, there is no need to assume asymptotic flatness in 
spacelike directions. This is important, because a self-gravitating thermodynamical 
system generically does not satisfy asymptotic flatness. In particular, the system 
semiclassically approximated by a black hole in equilibrium with Hawking radiation 
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is not asymptotically flat when the back-reaction of radiation on the geometry is 
taken into account. [5] Another advantage of using finite spatial boundaries appears 
in the treatment of rotation. Since any system in thermal equilibrium must rotate 
rigidly (if at all), [6, 7] such systems necessarily have finite spatial extent. Thus, a 
Kerr black hole surrounded by Hawking radiation can be treated as the semiclassical 
approximation to a thermal equilibrium system only if a spatially finite boundary is 
employed. As a final example, observe that the usual thermodynamic limit requiring 
infinite spatial extent does not exist for an equilibrium self-gravitating system at 
nonzero temperature. This is because the system is unstable to gravitational collapse, 
or recollapse if a black hole is already present. The instability of such a spatially 
infinite system at fixed temperature is reflected in a formally negative value for the 
heat capacity, which, in turn, implies that the canonical partition function diverges. 
(See, for example, Ref. [8].) On the other hand a spatially flnite system can avoid 
such difficulties. For the gravitational fleld at relatively low temperature, the system 
is approximated semiclassically by flat space flUed with dilute gravitational radiation; 
at relatively high temperature the system is approximated semiclassically by a large 
black hole surrounded by sparse gravitational radiation. [5, 9] 

Self-gravitating systems in thermal equilibrium are typically spatially inhomo- 
geneous because of gravitational "clumping". In particular, the temperature of an 
equilibrium system may vary in space due to gravitational redshifting.[10] As a con- 
sequence, such systems are characterized not by a single temperature, but can be 
described by a temperature field which is a local function defined on the spatial two- 
boundary. [11, 12] Correspondingly, the thermodynamical conjugate of inverse tem- 
perature is not simply the total energy, but rather an energy surface-density which 
is a local function on the spatial two-boundary. The microcanonical or canonical 
descriptions of a self-gravitating system are obtained by fixing the energy surface- 
density or surface temperature (respectively) as boundary data. Generally, all ther- 
modynamical intensive and extensive variables are functions defined on the spatial 
boundary. The appropriate definitions of energy density as well as angular momen- 
tum density are discussed in detail in Ref. [13], and are reviewed in Sec. 3. 

The density of states for the gravitational field is defined here as a functional 
of the energy surface-density, momentum surface-density, and the two-metric on 
the spatial boundary of the system. It is expressed formally as a functional integral 
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over Lorentzian metrics satisfying the boundary conditions, and includes contribu- 
tions from manifolds of various topologies. We evaluate the density of states in a 
"zero-order" approximation in which the functional integral is approximated by its 
integrand evaluated at an appropriate saddle point. When the boundary conditions 
are chosen such that the system is approximated classically by a stationary axisym- 
metric black hole, then in the zero-order approximation the entropy (identified as 
the logarithm of the density of states in absolute units G = c= h = kB = l) equals 
1 /4 the area of the black hole's event horizon. This result applies to any stationary 
axisymmetric black hole, including those that are distorted relative to the standard 
Kerr family by stationary external matter fields. This result also extends to black 
holes coupled to electromagnetic and Yang-Mills fields. [4] (The result also does not 
appear to depend on axisymmetry.) When the boundary conditions for the density of 
states are chosen such that the system is approximated classically by flat spacetime, 
it is shown that the entropy vanishes in the zero-order approximation, as expected. 

In nonrelativistic mechanics, the canonical partition function is defined by a 
sum over energy levels weighted by the Boltzmann factor and appropriate degener- 
acy factors. In the cases we shall treat, this is generalized and expressed as a (func- 
tional) integral transform of the density of states. At the level of thermodynamics, 
the change of boundary data amounts to a (functional) Legendre transformation be- 
tween the energy density and the inverse temperature, which are thermodynamically 
conjugate variables. At the level of dynamics, this change of boundary data amounts 
to a canonical transformation and the energy density and inverse temperature are 
given by the boundary values of a canonically conjugate pair of variables. For this 
intepretation, canonical conjugacy is defined with respect to the history of the spatial 
boundary, not with respect to the usual spatial time slices. Analogous relationships 
hold for the angular momentum density and its conjugate, the angular velocity, as 
well as other pairs of conjugate variables. These results reveal an intimate connection 
between thermodynamical and canonical conjugacy for self-gravitating systems. [14] 

In Sec. 2, we present a real-time functional integral expression for the density 
of states in nonrelativistic mechanics. The relevant action functional is Jacobi's ac- 
tion, [15, 16] in which the energy of the system is fixed. Details of the construction 
are given in an Appendix. In Sec. 3, we draw on the analysis of Ref. [13] to obtain 
a "micro canonical action", an action functional for which the appropriate bound- 
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ary conditions include fixed energy surface-density, momentum surface-density, and 
boundary two-metric. The microcanonical action is used in Sec. 4 to express the 
density of states formally as a functional integral. The functional integral is then 
evaluated in the saddle point or zero-order approximation to show that the entropy 
of any stationary axisymmetric black hole is 1/4 the area of its event horizon. In 
Sec. 5, the canonical and grand canonical partition functions are derived from the 
microcanonical construction and the correspondence between thermodynamical and 
canonical conjugacy is described. The first law of thermodynamics is derived in Sec. 6 
by considering variations of the microcanonical action with respect to the boundary 
data. 

II. DENSITY OF STATES IN NONRELATIVISTIC 

QUANTUM MECHANICS 

In this section, the formal expression of the density of states as a real-time func- 
tional integral is derived for nonrelativistic systems with a finite number of degrees 
of freedom. Our starting point is the density of states expressed as 

u{E) = T:t5(E-H) , (2.1) 

where H is the Hamiltonian operator for the system. The number of quantum states 
between Ei and E2 is 

[ ^ dEv{E) , (2.2) 

J El 

as seen by taking the trace in a basis of energy eigenstates. Using a coordinate basis 
for the trace, the density of states becomes 

v{E) = J dx {x\5{E - H)\x) , (2.3) 

where x represents a set of configuration coordinates, a;^, a;^, The matrix ele- 
ments in the integrand above are the diagonal entries of the matrix {x"\5{E — H)\x')^ 
which can be expressed as 

1 r+00 

{x"\5{E - H)\x') = / dTe'^^/^{x"\e-'''''/^\x') . (2.4) 
2nh 
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In turn, the matrix elements in the integrand of Eq. (2.4) can be expressed as a 
functional integral, [1] 



~x(T)=x" 

{xy-'^'^/^lx') = / VHe'^^/^ . (2.5) 

Jx(0)=x' 



rx{T)=x" 
'x{0)=x' 

This functional integral is a sum over histories x{t) that begin at a;(0) = x' and end 
at x{T) = x", with VH denoting a measure for the space of histories. The histories 
are weighted by exp{iST/fi), where St[x] is the usual Hamilton's action with fixed 
time interval T. 

Collecting together the above results, the density of states becomes 

^^E) = / dT dx p^g.(s.+£;T)/ft _ ^2.6) 

27rft J Jx{0)=x 

This expression for z/(i?) is an integral over all histories x{t) that are periodic for 
some real time interval. In the Appendix, it is shown that this functional integral is 
precisely the sum over periodic histories constructed from Jacobi's action 5"^. [15, 16] 
In Jacobi's action, the energy E is fixed rather than the time interval T. Consequently 
the path integral for the density of states can be written as 

v{E) = / VHpe'^^^^ , (2.7) 



where it is understood that the histories Hp contributing to this path integral are 
periodic in real time. This is the key result that will be generalized to the case of self- 
gravitating systems: the density of states is a sum over periodic histories, weighted 
by a phase that is given by the action appropriate for describing the system at fixed 
energy. 

The canonical partition function is obtained by summing the Boltzmann factor 
over each energy level. In terms of the density of states, the partition function is 
given by a Laplace transform, 

/•oo 

Z{(3) = / dEy{E)e-f^^ , (2.8) 

JO 

where I3~^ = fcs x (temperature) and is Boltzmann's constant. Using the expres- 
sion (2.1) for i'{E) and assuming the energy spectrum is positive gives the familiar 
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result Z{P) = Tr ex.p{—HP). Alternatively, with the density of states expressed as 
the path integral (2.6), the partition function becomes 

Z{(3)= dEe-f^^-^ dTe'^^/^ VHpe'^^/^, (2.9) 
Jo 27r?i J_oo J 

where again Hp refers to periodic histories. With the change of variables T = —it, 
the partition function becomes 



rOO -1 ptOO / p 

Z{(3) = / dEe-^^ / d{T/n)e^^/^[ VEpc'^^/'' 

Jo 27^^ \J 



(2.10) 



This expression is simply the Laplace transform of the inverse Laplace transform of 
the functional integral in parenthesis. To be precise, the identification of the integral 
over r/h with an inverse Laplace transform assumes that the integration contour 
passes to the right of any poles in the complex plane. Such points will be ignored 
in the present formal analysis. Then the result of the successive inverse Laplace and 
Laplace transforms is to set r equal to in the path integral factor, leaving 



(2.11) 

-ihl3 



This is Feynman's result, [1] that the canonical partition function can be written as 

an "imaginary-time" path integral. 

The expression (2.11) for Z{P) is often taken as the starting point for a treat- 
ment of thermodynamics by functional integral methods. Observe that the existence 
of the canonical partition function depends on the convergence of the Laplace trans- 
form (2.8). If the density of states increases too rapidly for large E, then Z{P) is 
not defined. This occurs when the heat capacity for the system is formally negative, 
signaling a thermodynamical instability. (The relationship between the convergence 
of the Laplace transform for Z((3) and the sign of the heat capacity is spelled out in 
Ref. [8].) We regard the real-time functional integral (2.7) for the density of states 
as a more fundamental expression. In the following sections, this result is generalized 
to the case of self-gravitating systems. 
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III. MICROCANONICAL ACTION 



We begin by summarizing the results of Ref. [13]. Start with the action for 
gravity 

t" 

S[g]^^f d^xyJ^{^-2K) + - I d^x\/hK--l d^x^/^e-S^, (3.1) 

where k = Stt and Newton's constant is set to unity. The spacetime manifold is 
M = S X /, the product of a space manifold E and a real line interval /. The two- 
boundary of space E is 5, and the history of B is^B = B x I. The submanifolds 

of M that coincide with the cndpoints of the line interval / are the hypersurfaces 

t" 

t' and t" . The notation J^, d^x represents an integral over t" minus an integral 
over t' . We also use the following notational conventions. The metric and curvature 
tensor on spacetime M are g^i, and 3?^z/a/3 respectively, the metric and extrinsic 
curvature on the hypersurfaces S are hij and Kij respectively, and the metric and 
extrinsic curvature on ^B are 7y and respectively. (Latin letters z, j, etc. are 
used as indices for tensors on both ^B and E. No cause for confusion arises from this 
convention.) The term in equation (3.1) is a functional of the metric 

lij on ^B- 

however, it will be seen that such a term is unnecessary. 

The action S is written in canonical form by foliating M into spacelike hyper- 
surfaces S. Without loss of physical generality, we restrict these hypersurfaces to be 
orthogonal to the boundary element ^B. That is, on the boundary element "^S, the 
timelike unit normal of each surface S is required to be orthogonal to the spacelike 
unit normal of ^B. The result is [13] 

S= I d'^x[P'^hij - Nn-V'Hi] - [ d^xV^[Ns-V'ji\ , (3.2) 

where N is the lapse function, is the shift vector, and the gravitational momentum 
conjugate to the metric hij is 

P'' = -^^{h^'h^' - h''h^'){hki - 2D,kVi,) . (3.3) 
4k iv 

The gravitational contributions to the Hamiltonian and momentum constraints are 

n=^ [2P'^Pi, - [Plf] - |5(i? - 2A) , (3.4a) 
= -2D J Pi , (3.46) 
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where R and Dj are the curvature scalar and covariant derivative on E, respectively. 
In the surface term of the action (3.2), a denotes the determinant of the metric tensor 
on S, and the energy surface-density e and momentum surface-density ji are defined 

by 

= (3.56) 

Here, Uy, rii and k denote (respectively) the induced metric, the unit normal, and 
the trace of the extrinsic curvature for S as a surface embedded in S. In writing the 
Hamiltonian action, we have assumed that 5°, if present, is a linear functional of the 
lapse and shift on ^S, in accordance with the discussion of Ref. [13]. 
The Hamiltonian obtained from the action (3.2) is 

H= I (fx[Nn + V'ni]+ [ d^x^/^[Ns-V'ji] . (3.6) 

The shift vector at the boundary B must satisfy njF*|^ = 0, so that the Hamilto- 
nian does not generate spatial diffeomorphisms that map the field variables across the 
boundary B of the space manifold S. This restriction implies that the Hamiltonian 
evolves the initial data into a spacetime whose foliation by spacelike slices is orthog- 
onal to the boundary element ^B. With the surface terms that appear in Eq. (3.6), 
the Hamiltonian has well defined functional derivatives with respect to the canonical 
variables under the conditions that AT, F*, and aab are fixed on the boundary B. 

We use indices a, b, etc. to denote components of tensors on B. Such tensors 
also can be viewed as tensors on S that are orthogonal to the unit normal rf of 
B. Thus, for example, we write the two-metric on B as aab or (7^, the extrinsic 
curvature of B embedded in E as kab or kij, the shift vector on B as V"' or V^, 
and the momentum surface-density as ja or jV We also have occasion to view these 
tensors as tensors on spacetime, and will then use spacetime indices n, etc. 

One can calculate, as in Ref. [13], that a general variation of the action S with 
respect to the canonical variables hij, lapse N, and shift is given by 

t" 

5S = (terms giving the equations of motion) + / d^xP^^Shij 

Jt' 

- j^^ Sx^ [e - jaSV^ - iN/2)s^^5aab] • (3.7) 



The term s**^ is the surface stress tensor on B, defined by[13] 



S 



[A;"'' + {ma' - - ^ , (3.8) 



where a* is the acceleration of the timehke unit normal of the spacelike hypersur- 
faces. The expression (3.7) shows that suitable boundary conditions for S are found 
by fixing the induced metric on the boundary dM. That is, fix the three-metric 
components hij on t' and t" , and fix the three-metric components N, V"', and aab on 
^B. Then the surface terms in the variation 5S vanish, and solutions of the equations 
of motion extremize the action S with respect to variations that obey these boundary 
conditions. 

What we define as the microcanonical action Sm is obtained from 5" by adding 
boundary terms that change the appropriate boundary conditions on ^B from fixed 
metric components A^, V, and cXab to fixed energy surface-density e, momentum 
surface-density ja, and boundary metric aab- Thus, define 

Sm = S+ f d^x^/^[Ne-V''ja] (3.9a) 
= I d^x[P'^h,j - NH-V'Ui] , (3.96) 

JM 

and from Eq. (3.7), the variation of Sm is 

t" 

5Sm — (terms giving the equations of motion) + / d^xP^-'Shij 

Jt' 

+ [ d^x\N6{V^e)-V''6{V^ja) + {NV^/2)s''^Saab] . (3.10) 

This result shows that solutions to the equations of motion extremize Sm under 
variations in which e, ja, and aab are held fixed on the boundary B. Observe that 
the unspecified subtraction term 5° does not appear in the action Sm., so in this 
sense the microcanonical action is unique. Nevertheless, the variation (3.10) of Sm is 
expressed in terms of the surface stress-energy-momentum components e, ja, and s"^, 
which do depend on S^ for their definitions. However, since S^ is a linear functional 
of the lapse and shift, the 5° dependences contained in the various terms of SSm 
actually cancel. 
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The boundary terms in the variation (3.10) of Sm show that N and ^/as are 
canonically conjugate, where canonical conjugacy is defined with respect to the 
boundary element ^B. Likewise, and —y/aja are canonically conjugate, as are 
{N^/a/2)s°'^ and aah- The boundary terms added to S in Eq. (3.9a) to obtain the 
microcanonical action Sm amount to the addition of terms of the form "pq'" at ^B. 
These terms have the effect of changing the appropriate boundary conditions from 
fixed "q'" to fixed conjugate "p" . 

The microcanonical action (3.9) can be written in spacetime covariant form by 
using expression (3.1) for S and the decomposition of the scalar curvature 

3? = i? + K^.K^^" - {Kf - 2V^iKu^' + a'^) . (3.11) 

The extra boundary terms in Sm are written covariantly by using the decomposition 
of the extrinsic curvature found in Ref. [13]. That analysis yields the relation- 
ships 

k = {g^'' + u^u'')Q^, , (3.12a) 
-2ViP'^nj/Vh = -V^u^Qi^^/k , (3.126) 

for the corresponding terms in e and ji (see Eqs. (3.5)). The microcanonical action 
in spacetime covariant form is therefore 

t" 

Sm[g] = — I d^x^m - 2A) + - / d^x^/hK -- I d^x^/^ta&'"'da . 
2k Jm « Jt' « Jm 

(3.13) 

Here, t is the scalar field defined on ^B that labels the foliation on which e, ja, and aab 
are fixed, 0'*^ is the extrinsic curvature tensor of ^B, and is the time vector field 
defined on ^B that specifies the time direction. In terms of the timelike unit normal 
u^^ of the slices B C^B, these quantities are given by = —Nd^t = {t^ — V^)/N. 

IV. MICROCANONICAL FUNCTIONAL INTEGRAL 

In Sec. II we showed that for nonrelativistic mechanics the density of states is 
given by a sum over periodic, real time histories, where each history contributes a 
phase determined by the action that describes the system at fixed energy. In the 
case of nonrelativistic mechanics, the energy is just the value of the Hamiltonian 
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that generates unit time translations. For a self-gravitating system, the Hamiltonian 
has a "many-fingered" character: space can be pushed into the future in a variety 
of ways, governed by different choices of lapse function N and shift vector F*. The 
value of the Hamiltonian (3.6) depends on this choice. More precisely, the value of the 
Hamiltonian is determined by the choice of lapse and shift on the boundary B, since 
the lapse and shift on the domain of E interior to B are Lagrange multipliers for 
the (vanishing) Hamiltonian and momentum constraints. Accordingly, the energy 
surface-density e and momentum surface-density ja for a self-gravitating system 
play a role that is analogous to energy for a nonrelativistic mechanical system. In 
particular, the energy surface-density s is the value (per unit boundary area) of the 
Hamiltonian that generates unit magnitude time translations of the boundary B, in 
the spacetime direction orthogonal to E. Likewise, the momentum surface-density 
ja is the value (per unit boundary area) of the Hamiltonian that generates spatial 
diffeomorphisms in the d/dx"" direction on the boundary B. 

The above considerations lead us to propose that the density of states for a 
spatially finite, self-gravitating system is a functional of the energy surface-density 
e and momentum surface-density ja- In addition to these energy-like quantities, the 
density of states is also a functional of the metric (Tab on the boundary S, which 
specifies the size and shape of the system. In the absence of matter fields, these 
make up the complete set of variables and v[s.,jaiO'ab\ is interpreted as the density 
of quantum states of the gravitational field with energy density, momentum density, 
and boundary metric having the values e, ja, and aab- The action to be used in the 
functional integral representation of v is Sm-, which describes the gravitational field 
with fixed £, ja, and aab- Note that e, ja, and aab play the role of thermodynamical 
extensive variables. These variables are all constructed from the dynamical phase 
space variables {hij, P*-') for the system, where the phase space structure is defined 
using the foliation of M into spacelike hypersurfaces. (We expect this to be a defining 
feature of extensive variables for general systems of gravitational and matter fields.) 
On the other hand, the variables and /2)s°'^ are not constructed from 

phase space variables. However, these variables are canonically conjugate to \/a"£, 
— \foja and aab where canonical conjugacy is defined with respect to the boundary 
element ^B. In Sec. VI, the relations of N ^ V^, and s"* to the intensive variables 
thermodynamically conjugate to ^/ae, —\faja, and aab are given. 
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By analogy with the functional integral (2.7) for the density of states in non- 
relativistic mechanics, the density of states for the gravitational field is formally 
expressed as 

iy[e,j,a] = ^ / 1^Hexp{tS^) . (4.1) 

(Planck's constant has been set to unity.) The sum over M refers to a sum over 
manifolds of different topologies. The three-boundary for each M is required to have 
topology dM = B X S^. If B has two-sphere topology, then the sum over topologies 
includes M = (ball) x S^, with dM = d{hall) x = S'^ x S'^ . Another example is 
M = (disk) X S'^, with dM = ^(disk) x S'^ = S'^ x S'^. The action Sm that appears in 
Eq. (4.1) is the microcanonical action (3.13) of the previous section, but with the t' 
and t" terms dropped because the manifolds considered here have a single boundary 
component dM — ^B: 

Smig] [ d^'x^i^ _ 2A) - - / d^x^t^&^^da . (4.2) 

Jm JdM 

The functional integral (4.1) for u is a sum over Lorentzian metrics g^^. Note that 
the action (4.2) may require the addition of a term that depends on the topology of 
M, such as the Euler number. 

In the boundary conditions on dM — B x , the two-metric aab that is fixed 
on the hypersurfaces B is typically real and spacclike. Likewise, the energy density 
e is real, which requires the unit normal to dM to be spacelike. Therefore, the 
Lorentzian metrics on M must induce a Lorentzian metric on dM, where the timelike 
direction coincides with the periodically identified S^. Note, however, that there are 
no nondegenerate Lorentzian metrics on a manifold with topology M = (disk) x 5'^ 
that also induce such a Lorentzian metric on dM. This implies that the formal 
functional integral (4.1) for the density of states must include degenerate metrics. 
(For a discussion of the role of degenerate metrics in classical and quantum gravity, 
see Ref. [17].) 

Now consider the evaluation of the functional integral (4.1) for fixed boundary 
data £, ja, Cab that correspond to a stationary, axisymmetric black hole. That is, start 
with a real Lorentzian, stationary, axisymmetric, black hole solution of the Einstein 
equations, and let T = constant be stationary time slices that contain the closed 
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orbits of the axial Killing vector field. Next, choose a topologically spherical two- 
surface B that contains the orbits of the axial Killing vector field, and is contained 
in a T = constant hypersurface. From this surface B embedded in a T = constant 
slice, obtain the data s, ja, and aab- In the functional integral for i/fe, j, cr], fix this 
data on each t = constant slice of dM. Observe that, to the extent that the physical 
system can be approximated by a single classical configuration, that configuration 
will be the real stationary black hole that is used to induce the boundary data. 

The functional integral (4.1) can be evaluated semiclassically by searching for 
four-metrics g^i, that extremize Sm and satisfy the specified boundary conditions. 
Observe that the Lorentzian black hole geometry that was used to motivate the 
choice of boundary conditions is not an extremum of Sm, because it has the topology 
[Wheeler (spatial) wormhole] x [time] and cannot be placed on a manifold M with a 
single boundary x S^. However, there is a related complex four-metric that does 
extremize Sm, and is described as follows. Let the Lorentzian black hole be given by 

ds^ = -N^dT^ + hij {dx' + V'dT) {dx^ + V^dT) , (4.3) 

where N, V^, and hij are T-indcpendent functions of the spatial coordinates x"^. 
The horizon coincides with = 0. For convenience, choose spatial coordinates 
that are "co-rotating" with the horizon. [18, 12] Then the proper spatial velocity of 
the spatial coordinate system relative to observers at rest in the T = constant slices 
vanishes on the horizon, (V^/N) = 0, and the Killing vector field d/dT coincides with 
the null generator of the horizon. [18, 19] By assumption, the metric (4.3) satisfies 
the Einstein equations, which are analytic differential equations in T. Therefore the 
Einstein equations are satisfied by the above metric with T imaginary, or equivalently, 
with the replacement T —iT. This leads to the complex black hole metric 

ds^ = -{-iNfdT^ + hij{dx' - iV'dT){dx^ - iV^dT) , (4.4) 

where the coordinate T is real. 

The complex metric (4.4) satisfies the Einstein equations everywhere on a man- 
ifold with topology M = (disk) x S"^ , with the possible exception of the points AT = 
where the foliation T = constant degenerates. The locus of those points iV = is a 
two-surface called the "bolt". [20] Near the bolt, the metric becomes 

ds^ ^ N^dT^ + h^jdx'dx^ , (4.5) 
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and describes a Euclidean geometry. The sourceless Einstein equations are not sat- 
isfied at the bolt if this geometry has a conical singularity in the two-dimensional 
submanifold that contains the unit normals to the bolt for each of the T = constant 
hypersurfaces. However, there is no conical singularity if the circumferences of circles 
surrounding the bolt initially increase as 27r times proper radius. The circumference 
of such circles is given by PN, where P is the period in coordinate time T. Therefore 
the absence of conical singularities is insured if the condition 

P{n'd,N) = 2tt (4.6) 

holds at each point on the bolt, where n* is the unit normal to the bolt in one of the 
T = constant surfaces. Because the unit normal is proportional to diN at the bolt, 
condition (4.6) restricts the period in coordinate time T to be P = 2ti/kh-, where 
Kh = [{diN)h'^^ {djN)Y/'^\_^ is the surface gravity of the Lorentzian black hole (4.3) 
(not to be confused with the constant k = Stt that appears in the action (3.1)). Note 
that the surface gravity of a stationary axisymmetric black hole is a constant on its 
horizon, [19] so the period P — 2TxjKn satisfies the condition (4.6) at each point on 
the bolt. 

The lapse function and shift vector for the metric (4.4) are = —iN and 
= —iV^. Thus, the complex metric (4.4) and the Lorentzian metric (4.3) differ 
only by a factor of —i in their lapse functions and shift vectors. In particular, the 
three-metric hij and its conjugate momentum P^^ (see Eq. (3.3)) coincide for the 
stationary metrics (4.3) and (4.4). [12] Since the boundary data e, ja, and aab are 
constructed from the canonical variables only, the complex metric (4.4) satisfies the 
boundary conditions imposed on the functional integral for i/[£,_7, cr]. 

The complex metric (4.4) with the periodic identification given by Eq. (4.6) ex- 
tremizes the action Sm and satisfies the chosen boundary conditions for the density 
of states z^[£,j, cr]. Although this metric is not included in the sum over Lorentzian 
geometries (4.1), it can be used for a steepest descents approximation to the func- 
tional integral by distorting the contours of integration for the lapse N and shift 
in the complex plane. Then the density of states is approximated by 

i^[e,j, a] « exp{iSm[-iN, -iV, h]) , (4.7) 
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where Smi—iN^—iV^h] is the microcanonical action (4.2) evaluated at the complex 
extremum (4.4). The density of states can be expressed approximately as 

i^[s,j,a] ^ exp{S[e,j,a]) , (4.8) 

where S[e,j,a] is the entropy of the system. Then the result (4.7) shows that the 
entropy is 

S[s,j,a] « iSm[-iN, -iV,h] (4.9) 

for the gravitational field with microcanonical boundary conditions. 

In order to evaluate Sm for the metric (4.4), first perform a canonical decompo- 
sition for the action (4.2) under the assumption that the manifold M has the topology 
of a punctured disk xS'^. That is, the spacelike hypersurfaces E have topology I x S'^, 
and the timelike direction is periodically identified (5*^). The outer boundary of the 
disk corresponds to the three-boundary element of M (denoted dM in Eq. (4.2)) 
on which the boundary conditions £, ja, and aab are imposed. The inner boundary of 
the disk, the boundary of the puncture, appears as another boundary element for 
M. (No data are specified at "^H.) The canonical decomposition is largely a reversal 
of the steps that lead from the form (3.9) for Sm to expression (3.13), which applies 
when S has a single boundary B. In the present case, boundary terms appear at 
from the volume integral of the term V ^{Ku^ -\-a^) in 3?, and from an integration by 
parts on the term involving the shift vector and momentum constraint. The result 
is* 

Sm= [ d'^x[P'^hij-Nn-V'ni]+ f d^x^/a[n'{^iN)/K+2niVjP'^ /^/h] , (4.10) 

where the expression Oj = (diN)/N for the acceleration of the timelike unit normal 
has been used. 

Now evaluate the action Sm on the punctured disk x S'^ for the complex metric 
(4.4), and take the limit as the puncture disappears to obtain a manifold topology 
M — (disk) X S"^. In this limit, the smoothness of the complex geometry is assured 
by the regularity condition (4.6). Since the metric satisfies the Einstein equations, 
the Hamiltonian and momentum constraints in Eq. (4.10) vanish, and the terms 

* The boundary term at has been given in Ref. [12]. 
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P'^^hij also vanish by stationarity. Moreover, the second boundary term at is 
zero because the shift vector vanishes at the horizon. Thus, only the first of the 
boundary terms at survives. Evaluating this term for the complex metric (4.4), 
that is, for the lapse function N = —iN, and using the regularity condition (4.6), the 
microcanonical action becomes 

Sm[-iN,-iV,h] = -^J dT j d^x\^n'diN 
= Ah 

K 

= -\Ah ■ (4.11) 

Here, A„ is the area of the event horizon for the Lorentzian black hole (4.3). 

The result (4.11) for the microcanonical action evaluated at the extremum (4.4) 
leads to an approximation for the entropy (4.9), which is 

S[e,j,a]^-^A, . (4.12) 

The generality of the result (4.12) should be emphasized: The boundary data £, j, 
and (7 were chosen from a general stationary, axisymmetric black hole that solves the 
vacuum Einstein equations within a spatial region with boundary B. Outside the 
boundary S, the black hole spacetime need not be free of matter or be asymptotically 
flat. Thus, for example, the black hole can be distorted relative to the standard Kerr 
family. Furthermore, recall that the quantum-statistical system with this bound- 
ary data must be classically approximated by the physical black hole solution that 
matches that boundary data. The result (4.12) shows that the entropy of the sys- 
tem is approximately 1/4 the area of the event horizon of the physical black hole 
configuration that classically approximates the contents of the system. 

It also should be emphasized that the microcanonical action Sm is independent 
of the term in Eq. (3.1); thus the entropy is independent of as has been shown 
in the framework of the canonical partition function. [9] Moreover, by setting 5° = 
in the definitions (3.5a) and (3.5b), the boundary data can be taken to be £ = k/K, 
ji = ~2aijnkP^^/Vh, and aab- 
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The calculations above have been carried out in the "zero-order" or classi- 
cal approximation. Beyond this approximation, the density of states will acquire 
a contribution arising from integration over quadratic terms in the functional inte- 
gral. Correspondingly, the entropy will acquire corrections to the zero-order result 
Afj/A. Physically, the system can be viewed in the zero-order approximation as 
consisting of a "vacuum" black hole. The next order contribution to the functional 
integral is viewed as arising from thermal gravitons surrounding the black hole. It 
is known that any stationary, axisymmetric system in thermodynamical equilibrium 
must rotate rigidly. [6, 7] Therefore the average distribution of graviton radiation sur- 
rounding the black hole must rotate rigidly with an angular velocity equal to that of 
the black-hole horizon. As a consequence, an equilibrium thermodynamical system 
cannot have infinite spatial extent, because the graviton fiux would then exceed the 
speed of light beyond some speed-of-light surface surrounding the black hole. This 
conclusion is supported by the analysis of Frolov and Thorne, [7] who show that for a 
quantum field in the Hartle-Hawking vacuum state on a Kerr black hole background, 
the Hadamard function is singular on and outside the speed-of-light surface. The 
above observations indicate that the density of states calculation of this section is 
not valid if the two-boundary B used to generate the boundary data £, ja-, and Uab 
is too far from the (rotating) black hole. The difficulty should show itself in the 
calculation of the quadratic contribution to the functional integral for the density 
of states. One possibility is that for a too-large boundary S, the contour for the 
functional integral cannot be distorted from Lorentzian metrics to pass through the 
extremum (4.4) along a path of steepest descents, but only along a path of steepest 
ascents. In this case, there may be no (generally complex) classical solution that 
dominates the functional integral for the density of states. 

Finally, consider the steepest descents evaluation of the density of states (4.1) 
for boundary data e, ja, and aab that correspond to flat Lorentzian spacetime. That 
is, use a two-boundary S in a stationary time slice of flat spacetime to induce the 
boundary data, then fix this data on each t = constant slice of dM. In this case, 
the same flat spacetime that motivates the boundary conditions can be periodically 
identified and placed on a manifold with boundary topology B x S^. It therefore 
constitutes a saddle point for the functional integral for u. More precisely, continu- 
ously many saddle points are obtained since the periodic identification can be made 
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with any proper period. Since these saddle points all arise in a topological sector 
with M = T,x S^, the action Sm can be written in the form of Eq. (3.9). This shows 
that Sm vanishes at each of these saddle points, so the entropy (4.9) vanishes in this 
"zero-order" approximation. 

V. CANONICAL PARTITION FUNCTIONS 

The canonical partition function characterizes a system that is open to exchange 
of energy with its surroundings, and has fixed inverse temperature. In the case of 
self-gravitating systems, the inverse temperature P is fixed on the boundary B that 
separates the system from its surroundings. Recall that P is not, in general, constant 
on B. The partition function is defined by an integral over energy densities ^/o■s, 

Zc[P,j,cr] = J 1^[V^£]j^[£,j,o-]exp(^- j d'^xy/aeP^ , (5.1) 

where the exponential factor arises from a product of Boltzmann factors for each 
point of B. Using the approximate identification (4.8) of entropy S as the logarithm 
of the density of states, Zc becomes 

Zc[P, j, (t]^ j V[^e\ exp {s[e, j, ^] - ^ rf'x^e/?^ . (5.2) 

The partition function can be evaluated approximately by performing the integration 
over ^/ae in a steepest descents approximation. The stationary point in ^/ae is given 
by the solution e{P) of the equation 

JF^ = ^ ' (^-^^ 

which will be recognized as a generalized form of the usual relation between the 
entropy of a system and its thermodynamic temperature. The approximation for the 
canonical partition function becomes 

lnZ,[P,j,a]^S[s{p),j,a] - [ d\y/^Pe{p) , (5.4) 

JB 

which expresses the Massieu function InZg as a (functional) Lcgendre transform of 
the entropy S. The expectation value of energy density is defined by 

5 In Zr 



6P 

^ j V[y/as]iy{y/as)exp(^- j (fx^eP^ . (5.5) 
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This integral can be carried out in a steepest descents approximation, with the result 
{^/as) ^/ae{(3). 

By inserting expression (4.1) for the density of states into Eq. (5.1), the canon- 
ical partition function can be written as 



= V / VHexpliSm - i / d^x^/aNe 

J \ JdM 



J2 [ T>Heycp(iS, 

M 



J dtN\^ = -if3 



J dtN\^ = -il3 

(5.6) 



From the discussion of Sec. 3, it is clear that Sc is the action appropriate for boundary 
conditions consisting of fixed two-metric cTab, fixed momentum density ja, and fixed 
lapse N on dM. The functional integral (5.6) is a sum over Lorentzian metrics 
with these boundary conditions. Furthermore, the gauge invariant part of N on 
the boundary, namely the proper distance J dtN\^, is analytically continued to the 
imaginary value —ip. The distance / dtN\^ denotes the proper length of curves in 
the boundary dM = B x that are orthogonal to the slices B and begin and end 
on the same slice. If it is possible to rotate the contours of integration for the lapse 
function (at each point of M) to the imaginary axis, then the functional integral (5.6) 
for Z(. becomes a sum over Euclidean metrics with aab, jai ^ind j dtN\^ = P fixed 
on dM. This prescription for the functional integral representation of the canonical 
partition function generalizes the results of Gibbons and Hawking [3] to allow for 
a finite spatial boundary and the effects of rotation. Likewise, Eq. (5.10) below 
generalizes their results for the grand canonical partition function. 

The inverse temperature P that appears in the canonical partition function 
is the thermodynamic temperature of the system. It is measured by the so-called 
"zero-angular-momentum-observers" (ZAM0's)[21] at B, that is, by observers at 
rest on the spacelike slices B C^B, and whose four- velocities were earlier denoted by 
■u^. Likewise the "chemical potential" defined below is the angular velocity oj of the 
system at B as measured by these same ZAMO's. See Refs. [12, 22] for discussions 
of P and (jj as ZAMO-measured thermodynamical variables. 

For the grand canonical partition function Zg, the system is open to exchange of 
momentum as well as energy. In the self-gravitating case, assume as in the previous 
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section that the fixed boundary metric aab is axisymmetric, and let 0" denote the 
axial Killing vector field on B. The momentum density in the 0" direction is ^/Gja4>'^i 
and its thermodynamical conjugate is /?a; with a; denoting the chemical potential. 
Below, ijj is identified as the angular velocity of the system in the direction with 
respect to local proper time on B. (In Ref. [12] this proper-time angular velocity was 
denoted by a;.) The grand canonical partition function is defined by transforming 
both from fixed energy density \fae to fixed inverse temperature /? and from fixed 
angular momentum density \f^3a'¥^ to fixed /3a;: 

xexp(- / d^xy/al3{e - ujacp'')] ■ (5.7) 



\ Jb / 
(As defined here, Zg is still a functional of the component jj_ of ja in the direction 
orthogonal to 0". With axisymmetric boundary data, j± can be simply set equal to 
zero. Alternatively, Zg could be defined to include an integral transformation of jj^ 
to a zero value of its conjugate.) With the density of states approximated by the 
exponential of the entropy S, the stationary point for a steepest descents evaluation 
of Zg is given by the simultaneous solution of Eq. (5.3) and 

In the zero-order approximation, the grand partition function (5.7) becomes 



JB 



(5.9) 



where e and ja(j)°' are functions of (3 and uj that solve Eqs. (5.3) and (5.8). Equation 
(5.9) expresses the Massieu function IrvZg as a (functional) Legendre transformation 
of the entropy <S. The expectation values of ^/ae and ^/^ja4>°' defined by deriva- 
tives of InZg with respect to (3 and /5a;, respectively. The results are approximated 
by the solutions of Eqs. (5.3) and (5.8). 

Combining the functional integral expression (4.1) for the density of states with 
the definition (5.11) for the grand canonical partition function yields 

Zg [/3, Pu;,a] = J2j'^H exp (^iSm - Sx^^ {e - a;j«0«] 



M 



J dtN\^ = -il3 and J dtV't'\^ = -iPu> 
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Here, V'^ is the component of the shift vector in the 0" direction. For axisymmetric 
boundary data with jj_ = 0, the action Sg is precisely the action S discussed in Sec. 3 
for which the two-metric aab, lapse AT, and shift are fixed on the boundary dM. 
In the functional integral (5.10), the gauge invariant distance J dtN\^ is fixed to the 
value —iP, and J dtV'^\^ is fixed to the value —if3uj. The quantity J dtV'^\-^ gives 
the amount of "twist" in the periodic identification of the boundary dM = B x S^. 
More precisely, note that the curves on dM that are orthogonal to the slices B and 
begin and end on a single slice need not close. Then J dtV^\^ equals the proper 
distance separating the initial and final points of such a curve, as measured along 
a trajectory of the Killing vector field 0", where = ^/a^V^. If the contours 
of integration for the lapse N and shift are rotated to the imaginary axis of 
the complex plane, then the functional integral (5.10) becomes a sum over a set of 
complex metrics with dab, J dt = /3, and / dtV'^\^ = /3uj fixed on dM. 

Recall that the lapse function and shift vector V"' are canonically conjugate 
to energy density -/ae and momentum density —\foja-, respectively, where canonical 
conjugacy is defined with respect to the boundary dM. The functional integral 
expressions (5.6) and (5.10) for the canonical and grand canonical partition functions 
show that the canonical and thermodynamical conjugates of the extensive variables 
v^ff, and —\/aja are related by 



dtN 
dtV^ 



= -iP , (5.11a) 

B 

= -i(3u . (5.116) 



Furthermore, consider the partition function that is appropriate when the system is 
open to fluctuations in the two-boundary metric aab, and define an intensive vari- 
able {Ps/a/2)p°'^ that is thermodynamically conjugate to aab- Writing this partition 
function as a functional integral shows that 

J dt (iVv^/2)s"''|^ = -z(/3V^/2)p"'' , (5.12) 

where {N ^ /2)s°'^ is the canonical conjugate of aab and s"'^ is the spatial stress 
tensor (3.8). Relations (5.11) and (5.12) show that canonical and thermodynamical 
conjugacy are intimately connected. [14] Specifically, the thermodynamical conjugate 
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of an extensive variable equals i times the boundary value of the time integral of its 
canonical conjugate. These relations also hold when matter is minimally coupled to 
the gravitational field, [4] and can be generalized straightforwardly to cases of non- 
minimal coupling. 

Now consider choosing boundary data for the various partition functions such 
that the complex black hole solution (4.4) extremizes the corresponding action. In 
this case, the lapse and shift are given hy N = —iN and = —iV"'. Equation 
(5.11a) shows that the inverse temperature (3 equals the proper length of a curve 
orthogonal to the slices B in the boundary dM = B x oi the complex black hole 
(4.4). [11, 12] Prom Eq. (5.11b) the chemical potential is a; = V'I'/N. Thus, u is 
the proper angular velocity of the Lorentzian black hole (4.3) in the 0" direction, 
as measured by the ZAMO's.[ll, 12] Similarly, Eq. (5.12) shows that p"'' equals the 
spatial stress tensor s"^ for the Lorentzian black hole (4.3). 

VI. THE FIRST LAW 

The first law of thermodynamics expresses changes in the entropy of a system 
in terms of changes in the extensive variables. In the "zero-order" approximation, 
the first law follows from the general variation of the microcanonical action Sm- That 
variation includes terms that yield the classical equations of motion, plus boundary 
terms that arise from integrations by parts. Those boundary terms are just the ones 
displayed in Eq. (3.10) for the boundary ^B. Thus, the variation in Sm is 

SSm = (terms giving the equations of motion) 

+ / d^x\N SiV^e) -V'SiV^ja) + iNV^/2)s''''5aab\ • (6.1) 

JdM L 

If the variations are restricted to those described by complex black hole solutions of 
the form (4.4) for difi^erent choices of boundary data (extensive variables) e, ja, and 
(Tab, then the terms giving the equations of motion vanish and the variation becomes 

S{iSm) = JdtJ^ d'x [N5{y/^e) - V''5{y/^ja) + (iVV^/2)s«^5aa6] . (6.2) 

Using the identifications (5.11) and (5.12) for the complex black holes and the ap- 
proximation S ~ iSjn, this variation becomes 

5S[e,3,a]^5{A,/A) 

d^x [/3 SiV^s) - (3uj 5{V^jar) + PiV^p''''/2)5aab] • (6.3) 
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This is the first law of thermodynamics for the gravitational field in a spatially finite 
region. It is seen to have the form dS = ''^pdE — PojdJ + PpdV" , familiar from 
standard thermodynamical treatments of nongravitating systems. In fact, if the 
boundary data are chosen such that is a constant on B, then the first term in 
Eq. (5.15) becomes pdE, where E = J^d'^x^/ae is the total (quasilocal) energy of 
the system. [13] If the boundary data are chosen such that f3uj is a constant on S, 
then the second term in Eq. (5.15) becomes —(3LodJ, where J = d'^x^ja(f>°' is 
the total angular momentum in the (f)°' direction. [13] Likewise, if the boundary data 
are spherically symmetric, then the third term in Eq. (5.15) becomes PpdA, where p 
is the surface pressure and A is the surface area of -B.[9, 14] However, it should be 
emphasized that these simplifications hold simultaneously only when the formalism 
is restricted to static, spherically symmetric systems. In order to treat a system that 
is classically approximated by, say, a distorted Schwarzschild black hole or a rotating 
black hole, it is necessary to consider boundary data that are not constant functions 
on the boundary surface -B.[ll, 12] 

APPENDIX: PATH INTEGRAL FOR JACOBI'S ACTION 

Consider a system described by the phase space x^, pi, x^, p2, ■ ■ ■ , and let a 
denote a parameter along the phase space path that increases monotonically from a' 
at one endpoint to a" at the other endpoint. Suppressing the indices on x and p, 
Jacobi's action reads[16] 

Se= [ da{xp- Nn{x,p)} , (Al) 

where N is a Lagrange multiplier and H{x,p) = H{x,p) — E is a constraint that sets 
the Hamiltonian H{x,p) equal to E. When varied with x{a') = x' and x{a") = x" 
held fixed, this action yields Newton's equations of motion with the restriction that 
the energy take the value E. The Lagrange multiplier N has the interpretation as 
the lapse in physical time, 

dt = Nda. (A2) 
Note that Jacobi's action is invariant under the gauge transformation 

6x = [x, eU] , 6p^[p,en] , 6N ^ i , {A3) 
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with e(cr) vanishing at the endpoints u', a" . This transformation is just the canonical 
version of reparameterization invariance, which reflects the arbitrariness in the choice 
of a path parameter cr. 

We wiU now construct the sum over histories associated with Jacobi's action 
and show that its trace is precisely the density of states i^{E). The gauge redundancy 
will be handled using BRST methods. [23] Let tt denote the conjugate to AT, so the full 
set of constraints becomes tt = and 7i = 0. Introduce the ghost coordinate C and 
momentum V associated with the constraint 7i = 0, and the ghost coordinate —iV 
and momentum iC associated with the constraint tt = 0. The ghosts C, "P, V and C 
are all anticommuting. The original phase space variables, Lagrange multiplier and its 
conjugate, and ghost variables constitute an extended phase space with fundamental 
Poisson brackets 

\p,x\ = -l, {AA) 
[7r,iV] = -l, (A5) 
[P, C] = -1 , (A6) 
[C,V] = -l. {A7) 

The theory is rank zero [23] since the constraints tt, Ti, have vanishing Poisson brackets 
with one another. As a result the BRST generator is a simple sum of constraints 
multiplied by ghost coordinates: 

n = -mv + nc . (^8) 

BRST transformations are defined by [ ■ , Vie] , with e an anticommuting parameter. 
In the extended phase space, Jacobi's action becomes 



Se= [ da{xp + N7r + VC + CP + [i;,n]] , 

Ja' 



{A9) 



where is an anticommuting gauge fixing function on the extended phase space. 
From the nilpotency of the BRST generator, [O, O] = 0, the action (A9) is seen to 
be invariant under BRST transformations with C{a') = and C{a") — 0. The path 
integral associated with Jacobi's action is now written as 

Z4x",x') = J VHe'^^'^ , (yllO) 
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where the measure T>H is the product over time of the Liouville measure on the 
extended phase space. The conditions 

x{a') — x' , x{a") — x" , 

7r(c7') = , 7r(c7") = , 

C{a') = , C{a") = , 

C{(t') = , C{a") = , (All) 

are BRST invariant and imply O(a') = = 0((t"), and thus constitute a consistent 
set of boundary conditions [23] that we will adopt for the path integral (AlO). 

The Pradkin-Vilkovisky theorem[23] states that the path integral Ze{x",x') is 
independent of the choice of gauge fixing function ^/J. For the purpose of evaluation, 
a convenient choice is ip = PN, so the path integral (AlO) becomes 

{x" ,x') = J VxVp VNVtt VCVV VCW 

X expj da [xp + Ntt + VC + CP - iPV - NH] | . {A12) 

With this choice of ifj the ghost contribution to the path integral decouples, and can 
be independently evaluated using any among a variety of techniques. One method 
is to recognize that the ghost path integral equals the determinant of the operator 
d^/da^ acting in the space of functions that vanish at a' and a". This determinant 
can be regularized, [24] yielding the result [a" — a'). The integration over tt in the 
path integral (A12) gives a formal infinite product (over a) of delta functions of N, 
restricting the lapse function A?^ to be a constant. Thus, the result of the VnVN 
integration is to leave a single integral dNo/27r?i over the constant value A^o of the 
lapse. 

Collecting together the above results, the path integral becomes 
Z^x", x') = ^'^ ^ J dNo J VxVp expj ^ J da [xp - NqH] | . {A13) 

Using the identification (A2), the argument of the exponent can be expressed as an 
integral over t, while the integration variable No{a" — a') is seen to equal the total 
time interval T = J dt. This leads to 

Z^x", ^') = ^J dTe'^^/^ j VxVp expj dt [{dx/dt)p - i?] | , {AlA) 
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where the definition Ti = H — E has been used. The functional integral over x and 
p contained in Eq. (A14) gives the matrix elements of the evolution operator (2.5), 
so comparison with Eq. (2.4) shows that the path integral for Jacobi's action is 

Z^{x",x') = {x"\6{E - H)\x') . (A15) 

The trace of this path integral is 

i^{E) = J dxZE{x,x) , {A16) 

the density of states. 

The above analysis shows that the density of states is the sum over periodic 
histories constructed from Jacobi's action. Observe that this identification assumes 
the range of integration for T is over all real values. Thus, the path integral Ze{x" , x') 
differs from the causal Green function,* which is obtained from expression (A14) by 
integrating over just positive values of T. [25] 

By integrating the total time T over all real values, the path integral for ^{E) 
consists of a sum over pairs of histories, where the members of each pair are weighted 
with opposite phases. To see this, consider a typical periodic history x{t), p{t), 
with period T > that contributes with phase ex.p{iS/h) to the path integral for 
the density of states. Another history that contributes to I'iE) is x{t) = x{—t), 
p{t) = —p{—t) with period T = —T < 0. As t decreases from to T, the history £(t), 
p{t) passes through the same sequence of configurations as obtained from the original 
history for t ranging from to T. Observe that the history 5(t), p{t) is closely related 
to the time reversed history x{t) = x{T — t), p{t) = —p{T — t), which has period T 
and consists of the original sequence of configurations taken in reversed order (for 
increasing t). Now, if the system is time reversal invariant, then the actions for the 
original history and the time reversed history are the same: denoting the Lagrangian 
by L, 

S = S= f dtL{x{t),p{t)) 
Jo 



dtL{i{t),p{t)) , {A17) 

* In this context, the causal Green function is the Green function for the time- 
independent Schrodinger equation defined by the Fourier transform of the retarded 
Green function for the time-dependent Schrodinger equation. 
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where the last equahty follows from a change of dummy integration variables. But 
the phase in the path integral associated with the history x{t), p{t) is determined by 
S = dtL{x,p), so that S = —S. Therefore the history x{t), p{t) and the original 
history x{t), p{t) represent the same sequence of configurations but contribute with 
phases of opposite signs to the path integral for viE). Consequently, each such pair 
of histories contributes to z^(-E') with phase 2cos(S'/7i), confirming that the density 
of states is real. 
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